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Using the well known boson mapping, we relate the transverse magnetic susceptibility of a system of flux 
vortices in f+f dimensions to an appropriately defined conductivity of a one-dimensional boson system. The 
tilt response for a system free of disorder is calculated directly, and it is found that a subtle order of limits 
is required to avoid deceptive results. 
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1. Introduction 

A d+1 dimensional vortex system with correlated dis- 
order can be rppdcled as a d dimensional time dependent 
boson system.EI This mapping has proved extremely use- 
ful for the study of vortex behavior for systems with dis- 
order correlated in one direction.cl Such pinning may arise 
from artificially created ion tracks, screw dislocations, 
and/or twin or grain boundaries!! Unfortunately, diffi- 
culties in characterizing the disorder lead to uncertaintir^ 
in the interpretation of experimental measurements olI 
In contrast, in the 1+1 dimensional system of flux vor- 
tices in large Josephson junctions, disorder can be fab- 
ricated with great control.Q In particular, it is easy to 
create samples with disorder appropriately correlated in 
one direction^ This freedom i-has encouraged both the 
experimental! and theoreticalDTu study of these 1+1 di- 
mensional systems. Researchers have also found 1+1 di- 
mensional vortex models to be a useful simplification for 
understanding the physics of the more complicated 2+1 
dimensional systemsiiJ'Eil __ 

In this paper, expanding on previous worked we dis- 
cuss features of the boson mapping for 1+1 dimensional 
systems. The basic mapping is briefly reviewed in sec- 
tion 2. As might be expected, the response of the vortex 
system to a transverse magnetic field (the tilt response) 
is related to a conductivity of the boson system. As is 
well known, one must be very .careful about defining a 
conductivity in one dimension£3 The relation between 
the vortex tilt response and a boson conductivity is de- 
rived in section 3. We will see, however, that this result 
is extremely deceptive due to a subtle order of limits. To 
illustrate this problem, as well as the method of treating 
it properly, we explicitly calculate the tilt response for a 
system with no disorder in section 4. 

2. Boson Analogy 

We consider a system of vortices confined to the x — y 



plane and oriented along the y axis, with disorder cor- 
related in the y-direction. In this case, such the pinning 
potential can be written as LTXa; ) and the free energy Tj^f 
for a system of N vortices istrQ 
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where Xi{y) describes the path of the i vortex, H c \ is 
the lower critical field, h(x, y) = H x (x, y)/H y (x, y) is the 
local slope of the applied magnetic field, e is the vortex 
stiffness, and V(x) is the vortex-vortex interaction. 

If periodic boundary conditions are imposed on the 
system in the y direction, the partition function for the 
system of flux lines is given by 
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where y goes from to L. The sum over P is a sum over 
all possible permutations of boundary conditions, such 
that for all i, Xi(0) — Xj(L) for some j. As discussed in 
Rcf. |l|, in the limit of h{y) — 0, this partition function is 
identical to an imaginary time Feynman path integral for 
the grand canonical partition function of a collection of 
bosons in 1+1 dimensions. The y coordinate for the flux 
line system is mapped to the imaginary time coordinate 
of the boson system. The chemical potential of boson 
system is given by e(H/H c u — 1), and the boson mass 
is given by the fluxline stiffness e. The temperature k^T 
of the vortex is mapped to Plank's constant fi for the 
boson system, and the length L of the vortex system is (3U 
with fi the inverse temperature {X/ksTbose) of the bose 
system. The boson pair potential is given by the flux 
line pair interaction U(|x; — Xj\), and the y-correlated 
disorder potential U(Xi) for the vortex system also maps 
directly to the bose system. 
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3. Tilt Response 

The above described mapping is particularly useful in 
allowing one to make contact with the large body of 
knowledge of interacting boson systems. To make this 
comparison more fruitful, we now derive a relation be- 
tween the tilt response of the vortex array and the con- 
ductivity of the interacting bosons. As will be discussed 
in an explicit calculation in section 4, it is necessary to 
keep careful track of the order of limits defining the dc 
response to avoid deceptive results. 

To probe the linear response of the flux line system, a 
small field in the x direction can be applied. To avoid 
problems with-43ux lines piling up at the x boundaries 
of the system,EHl the field is applied over a finite width 
< x < W, and we will take the limit W — > oo at the 
end. For simplicity, we will choose H x to be given in the 
form 



h(x,y) = 



_ / h(y) for < x < W 







otherwise 
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although it is a trivial generalization to include more 
complicated variations of h(x, y) in the x-direction. The 
Fourier components of the function h(y) are defined by 
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where q must be of the form 2irn/L with n an integer. It 
should be noted that since h(y) is real, h(q) — h(—q)*. 
The tilt modulus t(q) can now be defined as 
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where 9(q; W) is the Fourier component of the local slope 



0(q; W) = \ I dye l ^e{y; W) 
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and 6(y; W) is the spatially averaged slope of the flux 
vortices in the range < x < W. More formally, 
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where s(x,y) is the local slope density 
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The expectation of the slope 6(q; W) is given by differ- 
entiating the log of the partition function 



(6(q;W)) = 



k B T d\nZ 
WLedh(-q)' 
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Thus, the tilt response is the second derivative, which 
can be written as 
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dye- iqy 6(0; W)6(y; W) ) (11) 



where we have used the fact that \ 9{q; W)j is zero in the 

limit of zero field, as well as the y-translational invariance 
of the system. In the limit L — > oo this expression be- 
comes analogous to the expression for tha-cpnductivity 
of the bose system at zero temperatureHLj and finite 
frequency w where the external field is applied over the 
range < x < W: 



a{u;W) 



itTiuj 



dTe -iur ( Tr j( T ) j( )) (12) 



where T T is the time ordering operator, and J(t) is the 
current at imaginary time r averaged over < x < W, 
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J ^ = W Jo dx ^ x,T ^ 
with the local current density 



dxi 



(13) 
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Thus with q taking the place of to, the formal relation is 

We ttH 



t(q;W) = qa(q;W) 
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4. Calculation of the Tilt Response: 
Clean Limit — Luttinger Liquid 

Eq.|l5| is appealing, but can lead to drastically incor- 
rect results if applied naively. In an interacting one- 
dimensional quantum system, we expect a finite dc con- 
ductivity. Taking the limit q — > 0, Eq.^| seems to im- 
ply a vanishing tilt response to a uniform field. It is 
intuitively clear, however, that the vortices should rotate 
uniformly to accommodate the transverse field if there 
are no pinning defects present. To resolve this apparent 
paradox, we now calculate the conductivity more care- 
fully for general q and W (with W much larger than the 
mean vortex/boson spacing but still finite). 

If there is no disorder in the system, the ground state is 
a periodic array of vortices. Fluctuations of this array are 
described by a dimensionless displacement field u(x,y), 
where x n (y) = £(n + u(n£, y)/2n), with £ the mean vortex 
separation. The free energy can then be written as 
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Here we have shifted the zero value of y. 

Within the boson mapping, can be interpreted as 

the bosonized action for the quantum particles. Generic 
values of K x and K y actually correspond to the non- 
Fermi-liquid "Luttinger liquid" state induced by interac- 
tions in one dimension. It is satisfying that this highly 
non-trivial behavior of the quantum system is encoded 
in our simple displacement field description. The values 
of the coefficients K x and K y are dependent on the pre- 
cise form of the interaction Vr-&s well as on the average 
spacing I. Standard methodsEj pan be used to evaluate 
these constants in various limits Ji-3 

By transforming into Fourier space, the free energy 
becomes a sum (ie, an integral) of uncoupled harmonic 
oscillator modes. The equipartition theorem then yields 
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where u qxtq is the Fourier transform of u(x, y). With 
the current (j(x, r)) or local slope density given in terms 
of the displacement field as s(x,y) = du/dy, the slope 
averaged from to W (see Eq.fi]) can be written as 
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Performing the x integration and substituting into Eq.|lT| 
yields the tilt response 
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By adding a small piece 8 2 to the (q x W) 2 term of the 
denominator, the integral can be performed by contour 
integration and the limit 8 — > can be taken at the end 
to give the result 



t(q;W) = 
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As noted above, this calculation can easily be generalized 
to the case where we are concerned with the response to 
a field h(x,y) that is an arbitrary function of x. 

Eq.^o] is an appealing physical result, and resolves the 
"paradox" described at the beginning of this section. To 
define the conductivity of the quantum system, we must 
take q — > for finite W, to prevent difficulties with equi- 
libration and transport across the ends of the sample. 
Then, using Imx^o f(x) = x/2 and Eq.[ll| the dc con- 
ductivity is finite. In the vortex system, however, it is 
clear that for finite W the fluxons cannot tilt, since they 
would have to pile up at the boundaries to match the un- 
tilted flux lines outside the field region. Taking therefore 



W — > oo first, the behavior f(x) — > 1 in this limit gives 
a finite t(q — 0). 



5. Conclusion 

We have shown that a surprising subtlety in the or- 
der of limits arises in the calculation of tilt response and 
conductivity even for the simple pure problem discussed 
above. The origin of the difficulty lies in boundary ef- 
fects, which are generally ignored in applications of the 
vortex-boson analogy. It is natural that transport prop- 
erties such as the conductivity are sensitive to changes 
at the boundary, due to the constraint of charge conser- 
vation. One might conjecture that in the localized Bose 
glass phase, the order of limits may become unimportant. 
This possibility, and the extent to which such problems 
manifest themselves in other types of correlation func- 
tions remain interesting open questions. 



Acknowledgments 
This research was supported by the National Science 
Foundation under Grant No. PHY89-04035 at the ITP, 
and National Science Foundation Grant DMR-9 1-1 5491 
at Harvard. 



1 D. R. Nelson, Phys. Rev. Lett. 60, 1973 (1988); D. R. 
Nelson and H. S. Seung, Phys. Rev. B 39, 9153 (1989). 

2 D. R. Nelson and V. M. Vinokur, Phys. Rev. B48, 13060 

(1993) . 

3 See, e.g. S. Flesher et. al., Phys. Rev. B47, 14448 (1993); 
M. Hawley et. al., Science 251, 1587 (1991); and Ref. f| 

4 R. H. Koch, V. Foglietti, W. J. Gallagher, G. Koren, A. 
Gupta, and M. P. A. Fisher, Phys. Rev. Lett. 63, 1511 
(1989); P. L. Gammel, L. F. Schneemener, and D. J. 
Bishop, Phys. Rev. Lett. 66, 953 (1991). 

5 L. Civale, A. D. Marwick, T. K. Worthington, M. A. Kirk, 
J. R. Thompson, L. Krusin-Elbaum, Y. Sum, J. R. Clem 
and F. Holtzberg, Phys. Rev. Lett. 67, 648 (1991); M. 
Leghissa, L. A. Gurevich, M. Kraus, G. Saemann-Ischenko, 
and L. Ya. Vinnikov, Phys. Rev. B 48, 1341 (1993). 

6 M. A. Itzler and M. Tinkham, Phys. Rev. B51, 435 (1994); 
unpublished; to appear in IEEE Trans. Appl. Supercond., 

(1994) . 

7 Leon Balents and Steven H. Simon, Phys. Rev. B, In Press. 

8 R. Fehrenbacher, V. B. Geshkenbein, and G. Blatter, Phys. 
Rev. B45, 5450 (1992). 

9 V. M. Vinokur and A. E. Koshelev, Zh. Eksp. Teor. Fiz. 
97, 976 (1990) [Sov. Phys. JETP 70, 547 (1990)]. 

10 T. Natterman, I. Lyuksyutov, and M. Schwartz, Europhys. 
Lett. 16, 295 (1991); J. Toner, Phys. Rev. Lett. 67, 2537 

(1991) ; Y.-C. Tsai and Y. Shapir, Phys. Rev. Lett. 69, 1773 

(1992) . 



3 



11 G. G. Batrouni and T. Hwa, Phys. Rev. Lett. 72, 4133 
(1994). 

12 D. S. Fisher and P. A. Lee, Phys. Rev. B 23, 6851 (1981). 

13 See, e.g. S. N. Coppersmith, D. S. Fisher, B. I. Halperin, P. 
A. Lee, and W. F. Brinkman, Phys. Rev. B 25, 349 (1982). 

14 For the physically relevant case of the Josephson junctions 
described in Ref. [?], near the lower critical field (H/ H c ur—i 
1 < 1) where the fluxon density is low, standard methodaij 
give (with notation from Ref. K x = ck&T '/ and 
K y — el/ (Auks Tc), with c a non-universal order one con- 
stant. In the high magnetic field limit (H 3> (j>o/dXj), these 
methods are not appropriate since the vortices overlap and 
become poorly defined. None the less, it can be shown 
that the free energy defined in Eq.^| still applies, with 
K x ^K y ^ ej/k B T. 

15 C. L. Kane and M. P. A. Fisher, Phys. Rev. Lett. 68, 1220 
(1992). 



4 



